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Abstract: Colour coherence in hard photoproduction is considered using the Monte 
Carlo event generators PYTHIA and HERWIG. Significant effects in the parton 
shower are found using multijet observables for direct and resolved photon induced 
reactions. The particle flow in the interjet region of direct processes shows a strong 
influence of string fragmentation effects. 

1 Introduction 

Colour coherence is an intrinsic property of QCD. Its observation is important in the study of 
strong interactions and in the search for deviations from the Standard Model ||]. It is interesting 
to look for colour coherence effects in hard photoproduction processes at the ep-collider HERA 
where large momentum transfers can be achieved and both direct (Fig. [l|(a)) and resolved 
(Fig. |l](b)) photon induced events occur. In Section multijet observables are studied which 




(a) (b) 



Figure 1: Examples of (a) direct and (b) resolved photoproduction 

reveal coherence at the parton shower level for both direct and resolved photoproduction. In 
Section |3| consideration of the particle flow in direct photoproduction shows colour coherence 
effects at the fragmentation stage of hadron production. 

2 Multijet observables 

The effects of colour coherence on the emission pattern of jets in e + e~ collisions are well known 
and intuitive. However in hadron-hadron collisions the large number of possible colour flows 
involved in jet production complicates the identification of variables sensitive to coherence. 
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Here, radiation patterns in 7p collisions are studied by considering events where soft radiation 
is hard enough to form a jet. This reduces the effect of secondary interactions in resolved 
photoproduction. 

The effects arising from different implementations of coherence were studied using 500 pb _1 
of events generated with PYTHIA 5.7 |2[ and HERWIG 5.8 ||. Direct and resolved events were 
generated separately and combined according to their cross sections. Events were generated 
with a minimum px of 20 GeV using the GRV proton and photon structure functions [|]]. Two 
event samples were generated using PYTHIA: the PYTHIA Coherent sample and the PYTHIA 
Incoherent sample which was obtained by switching off the coherence in the parton shower and 
the initial-final state coherence. HERWIG represents an alternate implementation of coherent 
processes. Jets of particles were found using the KTCLUS M algorithm in covariant pt mode 
with radius equal to 1. Three jet events with at least two jets of transverse energy (E\ et ) 



satisfying Ef 1 > 30 GeV and a third jet of Ef 1 > 3 GeV were selected. The jets are ordered by 
E{ et decreasing and referred to as "first" , "second" and "third" jet acccordingly in the following. 
Two scenarios were considered, one to reflect the acceptance in jet pseudorapidity (if 61 ) of 
the present ZEUS detector, \if et \ < 2.5, and one to show the possibilities with an extended 
acceptance, \if et \ < 4. In addition the events satisfied 0.2 < y < 0.85 and P 2 < 4 GeV 2 , where 
P 2 is the negative of the four-momentum squared of the photon. 

An overall drop in cross section is observed between incoherent and coherent event samples. 
For example, with a luminosity of 250 pb _1 and the standard detector acceptance, \rf et \ < 2.5, 
2600 multijet events are predicted by PYTHIA Incoherent, 1728 by PYTHIA Coherent and 
1665 events by HERWIG. For comparison in the extended acceptance scenario, \rf et \ < 4, 3012 
multijet events are predicted by HERWIG. 

The angular distribution of the third jet is also affected. Following || the angle a is defined 
as the azimuthal angle of the third jet about the second jet in the rj — if plane. Here, however, we 
use centre-of-mass (cm.) variables so a = arctan(Ai//| Aip\) where AH = sign (77^™) (773 — 1)2) 
and r]^ 71 = T] 2 — l/2(il2 — Vi) an d Aip = <p 3 — ip 2 . 771, 1)2 and i] 3 refer to the pseudorapidities in 
the lab frame of the first, second and third jets respectively and positive i] is in the direction 
of the incoming proton. A(p is the difference in azimuth ((f) between the second and third 
jets (in radians). The definition of a is illustrated for a typical event geometry in Fig. 0(c). 
The distribution of a as shown in Fig. |^(a) is broader for coherent events. This is consistent 
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Figure 2: a distributions (a) with the present acceptance and (b) with extended acceptance, (c) 
illustrates the definition of a for a typical event geometry. Jets are ordered by E\ et decreasing. 



with our understanding that for coherent processes radiation is generally suppressed in regions 
far from the directions of the incoming coloured partons. In addition, reducing the bias on 
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the distribution by increasing the acceptance from \r] 3 \ < 2.5 to \r] 3 \ < 4 produces a more 
pronounced depletion in the central region for coherent events (Fig. 0(b)). 

Canonical detector effects were simulated by smearing the HERWIG jet quantities with 
Gaussian distributions of varying widths. A resolution of 20% (10%) was used to smear the 
E\ et of jets with E\ et < 10 GeV (E{ et > 10 GeV). The width of the difference between generated 
and detected values of rf et and ip jet was taken to be 0.1. As shown in Fig. |2|(a) such detector 
effects should not seriously hinder the measurement of a distributions. 

The coherent emission of soft radiation does not have a strong effect on the jet profiles of 
the first and second jets. For instance, in Fig. [5](a) the transverse energy profile of the second 



jet is shown. This is the distribution of 5r/2 = rj 
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where n part is the r] of a particle within 
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Figure 3: Jet profiles for (a) the second and (b) the third highest E\ et jets for the extended 
acceptance scenario. The definition of the jet profile is illustrated in (c) for the second jet. 

one radian of ip of the jet centre, weighted by the transverse energy of the particle, as illustrated 
in Fig. 0(c). For this the extended acceptance scenario particles are considered with absolute r\ 
up to 5. The profile of the third jet is shown in Fig. 0(b). The occurrence of two peaks outside 
the jet core is due to partial overlap in ip of the first or second jet. A strong effect of coherence 
is apparent; it leads to less energy in the core of the third jet. 

One of the anticipated effects of colour coherence is that radiation from an incoming par- 
ton should be inhibited in regions far from the initial partons direction. Therefore in direct 
photoproduction events, where the single coloured parton in the initial state has positive rj, 
the coherent emission pattern should be at relatively higher rj than the incoherent emission. 
We have selected a subsample of events which is enriched in direct photon events by requir- 
ing x 1 > 0.8 where x 1 = (J^jets E 3 t e e~^ e ) / {2E 1 ) . The sum runs only over the two highest 
Ej? jets and E 1 is the energy of the incoming photon. The r\ of the third jet in the cm. 
frame, r]^ 71 = 773 — 1/2(772 — ?7i), is shown for this selection in Fig. |](a). The expected effective 
enhancement of radiation at large r)^ m can clearly be seen. 

Resolved events, with incoming coloured partons from both the 7 and p directions, should 
show an effective enhancement of radiation in coherent processes both at higher positive and 
at more negative pseudorapidities in comparison to incoherent emission. This effect is evident 
as shown in Fig. |](b). Note that the extended acceptance scenario must be employed in order 
to see the relative enhancement of radiation at high rj in resolved events. 

To summarize this section, a high integrated luminosity (~ 250pb _1 ) is desirable in order to 
accumulate statistics in multijet events at high E\ et . However luminosity upgrades which involve 
a significant reduction of forward acceptance are not worthwhile for this study. They destroy 
the sensitivity to colour coherence without significantly improving the statistical uncertainty. 
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Figure 4: r/g m distributions separated by an x 7 cut o/0.8 into (a) direct and (b) resolved samples 
for the extended acceptance scenario of \rf et \ < 4. In (c) the definition of r/f 11 is illustrated. 



3 Interjet string effects in direct photoprocesses 

Colour coherence effects which lead to a change in particle flow iV distributions in the inter- 
jet region should be rather pronounced in the direct photon induced processes such as QED 
Compton on quark (QEDC), QCD Compton (QCDC) and Photon Gluon Fusion (PGF). These 
distributions are considered here using the PYTHIA |2j generator with string (SF) or indepen- 
dent parton fragmentation (IF) into hadrons. Using SF is equivalent to taking into account the 
coherence effects at the hadronization phase of event generation. The flow iV depends on the 
string topology and colour antennae which are different for the three direct processes as shown 
in Fig. |(left). 

The calculation at the generator level was done using the GRV proton structure function 0] 
and minimum equal to 2.0 GeV. The HI detector simulation was taken into account as 
well. A jet-cone algorithm with radius equal to 1 was used for the selection of two jets 
and gamma-jet events with E 3 t etori > 3 GeV and jet (or final 7) emission angles 25 — 155°. 
This procedure corresponds to the selection of mainly direct processes. The calculated flow of 
charged particles with p t > 0.2 GeV emitted at angles of less than 20° to the reaction plane is 
shown in Fig. ^J(right) as a function of the scale angle Q. Q is defined || as the ratio of the 
particle angle 9h to the angle between partons. Q — corresponds to the direction of the initial 
state photon; Q = 1 - the final state quark; Q — —1 - the final state 7, gluon or antiquark for 
QEDC, QCDC or PGF respectively; Q = —2 or 2 - the proton remnant. 

It is seen that in the scale angle region between 1 and 2 (region II in Fig. [||) SF (solid 
histogram, generator level) and IF (dashed histogram, generator level) give different predictions 
for N. SF taking into account colour forces leads to a suppression of particle flow which 
is especially strong for QCDC process. The HI detector simulation (dark circles, SF) weakly 
distorts the generator level N distribution except for directions close to remnant proton emission 
where detector acceptance is rather low. Thus colour coherence effects can be observed at the 
detector level. 

It is interesting to consider ratios of particle flows N for different processes since the ra- 
tio is less sensitive to experimental errors. The ratios R = N(QCDC)/N(QEDC), R* = 
N(QCDC)/N(PGF) are shown in Fig. || Fig. ^ displays more clearly the role of colour co- 
herence which leads to drag effects in particle distributions. It is seen that for the case of SF 
the suppression in N for QCDC is more pronounced than for QEDC and PGF. At Q — 1.7 
the suppression reaches a factor of ~ 3. It has been found that misidentification of quark and 
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Figure 5: Left: Topology of direct photoproduction processes in the 7 — partem cm. frame. The 
double line is the proton remnant, dotted - strings. 9 q and 6h are quark and hadron emission 
angles. Right: Charged particle flow N normalized to 1 at the quark emission angle vs the 
scale angle Q. Solid line - generator level SF, dashed line - generator level IF, dark circle - HI 
simulation SF. In the regions I - IV Q is changing within the limits: < Q < 1 (I), 1 < Q < 2 
(II), -1 < n < (III), —2<Q<—1 (IV). 

gluon jets for QCDC does not change this conclusion. 

To observe colour coherence at the fragmentation stage of hadron production in direct 
processes it is necessary to distinguish these processes from resolved photoproduction and to 
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Figure 6: Ratios R and R* in the scale angle region 1 < Q < 2. Notations are as in Fig. [|. 

separate QEDC, QCDC and PGF from each other. The jet selection procedure used here 
enriches the data sample with direct processes. Further enrichment can be achieved by going 
to larger Ej etor ^ and by choosing events with x 7 close to 1. Since the direct processes cross 
section is low we expect ~ 230 QEDC events at the detector level for an integrated luminosity 
of ~ 100 pb -1 . So higher luminosity is needed to study interjet coherence. 

4 Conclusions 

The observation of colour coherence in photoproduction processes is an important challenge, 
particularly given the unique opportunity at HERA to study direct as well as resolved photon 
induced reactions. Since the cross sections of multijet events or of prompt photon reactions are 
small high luminosity ep-collisions are necessary for their study. 250 pb -1 would appear to be 
barely sufficient for these studies; however in upgrading to 1000 pb -1 it is essential that the 
forward acceptance should not be reduced. 
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